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Abstract 

We study the inclusive semileptonic rare decay b —* sl + l~ in minimal super- 
gravity model (mSUGRA). If tan/3 is large, down- type quark mass matrices 
and their Yukawa couplings cannot be diagonalized at the same basis. This 
induces the flavor violating neutral Higgs boson couplings. These couplings 
contribute significantly to decay b — ► sfi + fj,~ and b — > st + t~, but negli- 
gible to b — > se + e~ decay because of its negligible m e mass. The ratio 
R = B(b — > s/i + fi~)/B(b — * se + e~) can be very different from its corre- 
sponding value in the Standard Model. We find that part of parameter 
space can accommodate a large R value, and that maximum R value can 
be larger than 2. We also present our results in b — * st + t~ decay channel. 
Although it can not be detected now, it is potentially a new channel for the 
future observation of new physics. 

PACS numbers: 13.20.He; 14.80.Cp; 14.80.Ly 
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I. INTRODUCTION 



Rare semileptonic B decays provide an extremely helpful tool to search new physics 
beyond the Standard Model (SM). New physics contributions, which enter through one 
loop radiative corrections, may be observed whenever the SM contributions absent or 
suppressed. Thus, the measurement of B — > X s l + l~ has a very good chance to reveal new 
physics beyond Standard Model. 

One way in which new physics may reveal its presence is for there to be a deviation 
from the Standard Model prediction for B — > X s l + l~ . In the absence of a deviation in the 
cross section, it is still possible that new physics will be show itself in the details of the 
B — > X s l + l~ signal. The detailed study of this signal can also provide important clues as 
to the nature of the new physics. 

In this paper, we suggest that a significant discrepancy between B — > X s e + e~ and 
B — > X s fi + fi~ could arise in some instances of the two Higgs doublet model. Indeed such 
a discrepancy would strongly suggest new physics that couples to the mass of the lepton. 
In this paper we discuss how observable signals of this sort could arise in the popular 
framework of minimal super gravity. 

The SM predicts the b — > se + e~ and b — > sfi + fi~ branching ratio for M;+/~ > 2m M to 
be [1,2]: 

B{b -> se + e-) = 6.5 x 1(T 6 (1) 

B (b -> sn+fi,-) = 6.2 x 1(T 6 (2) 

B (b ^st + t-) =4.4 x 1(T 7 (3) 
So that the ratio of branching fractions, Rsm, is: 

B(b^ Sf i + fi-) 



Rsm = 



B(b^se+e-) 



~ .95 (4) 

SM 



SM predicts Rsm ~ 1 unless it is significantly altered by new physics. Thus the measure- 
ment of R provides a signal for specific classes of new physics. 

Supersymmetry (SUSY) is a promising candidate for new physics beyond Standard 
Model [3]. In SM, the decay b — > sl + l~ occurs through electroweak penguins and box 
diagrams. SUSY introduces several additional classes of contributions showing in Fig. 1: 

a gluino, down-type squark loop, 

b chargino, up-type squark loop, 

c chargino, up- type squark loop, (Higgs field attaching to charginos.) 

d neutralino down-type squark loop. 
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There are many papers exploring the photon penguin, Z penguin, gluino penguin and 
box diagram contributions in SUSY [4-6]. In our paper, we focus on the penguin contribu- 
tions to b — > sl + l~ mediated by neutral Higgs bosons. At the one-loop level, couplings of 
the "up-type" Higgs field H u to down-type quarks [7-9] are induced. This coupling gives 
a new contribution to the down type fermion mass matrix, and induces flavor violating 
couplings of neutral Higgs bosons. This induced flavor violating coupling increases with 
tan/3. If tan/3 is sufficiently large, this flavor violating coupling has significant contribu- 
tion to the branching fraction of B s — > fi + fi~ [10-20] and b — > sl + l~ decays [1,5,6,21-25]. 
In this paper, we would like to emphasize the contributions mediated by neutral Higgs 
exchange. As we shall see, the ratio R is particularly sensitive to new physics which 
contributes through this kind of mechanism, because other contributions which couple to 
e and /i equally contribute to both the numerator and denominator of Eq.4, therefore 
giving little contributions to R. 

Higgs mediated SUSY contributions to the amplitude increase with lepton mass mi. 
Therefore such contributions to b — > se + e~ are negligible while, given a sufficient large 
tan/3, the contribution to b — > sfi + fi~ may be comparable to the SM and so observable 
deviations of R from Rsm are possible. We will explore these effects as a function of 
SUSY parameter space which, in some places, are even larger than the SM contributions. 
In such cases, the ratio R = g^Ilge+e-) a ^ers substantially providing a possible way to 
detect new physics beyond SM. 

The main object of this paper is to examine the prospects for observing a deviation in 
R within the framework of minimal supergravity (mSUGRA) [26]. In Sec. II, we review 
the effective Hamiltonian for quark transition b — > sl + l~ . A detailed formula to calculate 
the branching fraction of b — > sl + l~ is presented. Sec. Ill contains a brief description of 
the mSUGRA model along with our main results. We analyze the ratio R in mSUGRA 
frame work and also discuss the large contribution of such effects to the decay b — > st + t~. 
We present our conclusions in the last section. 



The effective Hamiltonian for b — > s is derived by integrating out the heavy degrees of 
freedom at the electroweak scale or above. This Hamiltonian [27,13] can be written as: 



II. EFFECTIVE HAMILTONIAN 
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(5) 



o 2 
o 3 



q=u,d,s,c,b 



Only 07,9,10 and Qi,2 are obviously relevant to b —>■ sl + l directly. 
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04 = {8 alll P L bp) ]T (Qp-fPLQa), 

q=u,d,s,c,b 

05 = (Sa^PM E fo-fPRQfl), 

q=u,d,s,c,b 

6 = (s a ^P L bp) ]T &>fP R q a ), 

q=u,d,s,c,b 

7 = £^m b (s a a^P R b a )F^, 

08 = ^m h {s a T^a, v P R b )G a ^, 

e 2 

e 2 

010 = -^(Sa-fPMih^l), 

Qi = ^(s L b R )(ll), 

Q2 = ^(s L b R )(Hl), (6) 

a, (3 being color indexes, a labels the SU(3) generators, and Pl,r = (1 =F7s)/2. Cj are 
Wilson coefficients. 

0\ and 2 are current-current operators, 3 ...0 6 are QCD penguin operators. The 
contributions of these four quark operators to b — > s/ + /~ are proportional to the tree level 
matrix elements of operators 7 , 8 and 9 at one loop level. The chromomagnetic dipole 
operator, Og, gives no contributions to b — > sl + l~ . As a result, the (9i ...06 contributions 
can be absorbed by appropriately modifying the Wilson coefficients of operator O-j and 
Og, which originated in the Z° and 7 penguin diagrams with external l + l~ pairs. The 
neutral Higgs couplings SUSY contributions are mainly through Qi and Q 2 - The effective 
Hamiltonian for b — > sl + l~ is thus: 

H = -^v: s V tb {-2iC e 7 fS sa^{m b P R + rn s P L ) q -bhn + C e 9 ff s Ll ,b L h»l 

+CioS L 7 A1 M7' 1 75^ + C Ql s L b R ll + C Q2 s L b R l~f 5 l] (7) 

The Wilson coefficient Cg^ includes leading-order (LO) and next-to-leading order 
(NLO) logarithms, while and Ci enter only at the NLO level. In addition to the 
SM contributions to these coefficients, there are several classes of SUSY contributions. 
The photon penguin diagrams contribute to C7 . Three types of diagram: the photon 
penguin diagram; the Z penguin diagram; and the box diagram, contribute to Cg . The 
C10 is induced by the Z penguin diagram and the box diagram. 

In the calculation of Cj^, we find that Cj^ can be quite different from its SM value. 
There is a one-loop diagram with internal stop and chargino which gives a large contri- 
bution when tan/5 is large. This stop-Higgsino diagram is proportional to the product 
of the top and bottom Yukawa coupling constant, m t rrib/ \sinf3cosf3) , which grows with 
tan/?. There are no such terms in the calculation of Cg and C w . The corresponding 
stop-Higgsino diagram is proportional to the square of the top Yukawa coupling con- 
stants, m 2 t /sin 2 (3, which does not grow for large tan/5. As pointed out by Ref. [5], SUSY 
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contributions to Cg (to;,) and Cio(m&) are very small and alter these coefficients by < 
5% over the whole parameter space on Cg^ (m b ) and Ci (m 6 ). Because of this, SUSY 
contributions on Cg and C w can be ignored. 

We must consider the SUSY contribution to C7 . We write, 

C e 7 ff (M w ) = C* M {M W ) + Cf SY {M w ) 

where Cf M (M w ) is given in Ref. [27], and Cf USY (M w ) is taken from Ref. [6,28] with mass 
insertion approximation. For complete calculation, see Ref. [5]. At the m b scale: 

C eff = v W23 C eff {Mw) + 8 {r]U/23 _ ^2^^ + ( £ ^ )C72(Mm) 

6 i=l 

where 77 = ^^y, C^M^), C 2 (M W ), hi and Oj can be found in Ref. [27]. The expressions 
of Cg ff , and C 10 can be found in Refs [13,27]. Within SM, they are: 

C e 9 ff = C 9 + Y(s), C 9 = 4.138, Ci = -4.221 (8) 

where s = q 2 /ml, and g denotes the invariant momentum of the lepton pair. The 
expression for function Y(s) , coming from the one loop contributions of operators 0\ - 
Oq, can be found in Ref [13]. The Wilson coefficients Cq 1 and Cq 2 , can be found in Ref. 
[9]. In terms of the Wilson coefficients, the differential decay rate is [21,29]: 



dT GWml 



ds 128tt 5 



s s 



\(\<tf f \ 2 + \Cn\*) 



+ ^ s _ s) A (l, ^, ^) + (1 + S - - s - /*.)(! + ^)) 
V s s s / 

+ (|C7 9 | 2 - |C 10 |>(1 + ^ " s) + ||C^| 2 4(1 + ^) 

3 s s 

^2(1 + ^)(1 - /2 S ) 2 - S(l + 14/i s + /^) _ g2 (1 + ^ 

+2C 7 e// J Re(C 9 e// )(l + ^) ((1 - ^sf - S(l + p.) 
+\c 2 Ql {\ +fia -s)(s- Afn) + \c 2 Q2 s{l +fi s -s) 

+Ci Cq 2V /^(1 - fi a - s) 



(9) 



In Eq.9, jii = mf/mf, and ji s = m 2 /m 2 . The function X(x,y,z) is: 

A(x, y, 2) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz 

Integrating out s from 4m 2 /m 2 to (1 — m s /m b ) 2 , the branching ratios of decay b 
sl + l~ are easily obtained by: 



1 r (l-m s /m b ) 2 /dY\ 

B(b^ sZ+Z-) = i/ hrr Us (10) 
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III. RESULTS 



The minimal supergravity model [26] has been a popular model for SUSY phenomenol- 
ogy. It provides a well motivated realization of the Minimal Supersymmetric Standard 
Model (MSSM). In this model, SUSY is broken in a hidden sector. The fields in the 
hidden sector interact with usual particles and their superpartners only via gravity. In 
this way, SUSY breaking is communicated to the observable sector of Standard Model 
particles and their superpartners. 

The mSUGRA framework assumes that at the GUT scale (M GUT « 2 x 10 16 GeV), 
all scalar fields have a common SUSY breaking mass m , all gauginos have a mass mi, 
and also all soft SUSY breaking trilinear scalar couplings have a common value A . The 
resulting MSSM will have various soft SUSY breaking terms unified at M GUT . The soft 
SUSY breaking parameters are evolved from Mqut to weak scale using renormalization 
group equations (RGE). Minimization of the Higgs potential gives a relation between fi 2 , 
m|, B and tan/3, therefore fi 2 can be determined by m| [30] and the B parameter is 
traded in favor of tan/5, so that the model with radiative electrical symmetry breaking is 
completely specified by the SM parameters together with: 

mo, mi, Aq, tanfi, sign(fi) 

To evaluate the mass spectrum of the MSSM resulting from mSUGRA, we use the ISAS- 
UGRA program from the IS A JET package [31]. With these sparticle masses and mixing 
angles, we can get the Wilson coefficients Cq 1 and Cq 2 . Then the branching ratio of 
b — > sl + l~ can be computed using Eq. 9 and 10. 

We calculated the branching ratios of b — > se + e~ and b — > sfi + fi~. We find that if 
tan/3 is small, SM contributions dominate and SUSY contributions are very small. As 
tan/3 grows, SUSY contributions become large. O-j gives a sizable SUSY contributions 
with a little dependent on the lepton mass. The Or and Og interfere term also gives 
a comparable contributions. They together alter the branching fraction around 15% to 
20%. Qi and Q2 terms are induced by Higgs-mediated contributions. Their contributions 
to b — > sl + l~ branching fraction increase with tan 6 (3 and in some parameter space, they 
are even bigger than SM contributions. The corresponding Wilson coefficients Cqi and 
Cq2 are proportional to the lepton mass so that the branching ratio is proportional to 
lepton mass square. Therefore, Higgs-mediated SUSY contributions in decay b — > sfi + fx~ 
are about 10 4 bigger than that in b — > se + e~ decay and alter the ratio R = B(b — > 
sfi + B(b — > se + e~) deviating from its corresponding SM value. We require that the 
lepton pair mass is larger than 2 m M to reduce the dependences of the branching fraction 
on the lepton pair mass in small s range. 

In Fig. 2, we show the dependence of it! = B(b — > s/j, + iJ,~)/B(b — > se + e~) on tan/3 
and m for A = and A = —300 cases. We fixed mi/ 2 = 300 GeV in all the frames. 
Frame a), is R vs. tan/3 with fi < and m = 300 GeV. Frame b) is R vs. m with 
H < and tan/3 = 45. Frames c) and d) are for \x > case. Frame c) is R vs. tan/3 with 
m = 300 GeV, while frame d) has tan/3 = 53 with R vs m plot. The solid (dashed) line 
is for Aq = —300 (0) GeV, and the dotted line is the corresponding R value in SM. Values 
of tan/3 (mo) larger than the corresponding values denoted by circles on the curves are 
where m^ falls below its experimental bound 113 GeV, which gives out the most stringent 
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constraint. In frame b), m h > 113 GeV for all the regions in the graph along A = line. 
In frames c) and frame d), is always larger than 113 GeV. The main reason to show 
this figure is to understand the behavior of the ratio R. Several features need to be noted: 

• The R value is significantly larger for negative value of fi. As explained in ref. [9], 
changing the sign of \x changes the denominator of xfc in Eq.10 of ref. [9], so that 
a suppression for positive fi changes to enhancement for negative fi. 

• For fi < case, in low tan/3 regions, SM contributions dominate, so that all curves 
are close the SM value: R ~ 0.95. Same thing happens in m larger than 400 GeV 
range. 

• Again in fi < case, as tan/3 is larger than 44, there is a sharp rise for both A = 
and Aq = —300. This is because Higgs-mediated SUSY contributions are expected 
to behave as tan 6 (3/m\ [9] 2 . The solid and dashed curves indeed show this behavior 
as long as tan/5 is large. Clearly then, SUSY contributions are dominated by Higgs 
mediated penguin when tan/5 is large. 

• R deviates significantly from the SM value when tan/5 larger than 45 in fi < 
case. Unfortunately, in the mSUGRA scenario, these ranges are excluded by the 
experiment bound on the Higgs mass, rrih > 113 GeV [32]. 

• As shown in frame c) and frame d), in fi > case, R changes very slowly over the 
whole parameter space. Higgs-mediated SUSY contributions are very small [9]. The 
fact that solid line and dashed line are almost overlap and parallel to the dotted 
line tells us a constant R value exists. We even increase the value of A to around 
300 GeV. Still we get a R value, which is very close to its SM prediction. This 
excludes the possibility to find a signal of new physics in fi > case. 

Fig. 2 illustrates that typically large contributions to R occur in the fi < scenario 
with tan/3 > 40. For fi > 0, R is almost constant over the whole parameter space and very 
close to the SM prediction, so that no room has left for new physics. Although fi < is 
generally thought to be disfavored by the determination of the muon anomalous magnetic 
moment by the E821 experiment [33], a conservative estimate [34] of the theoretical error 
suggests that there is a region allowed [35] by this constraint, though perhaps in conflict 
with B(b — > S7) 3 , where b — > sl + l~ may provide the first hint of new physics if tan/5 > 40. 
This region would expand as more experimental data are accumulated. 



2 The contribution from h exchange is very small as long as h is a SM-like Higgs boson, and 
ran ~ ttla- 

3 It is worth reminding the reader that SUSY contribution, special for large values of tan/3, may 
have considerable theoretical uncertainty. Unlike constraints from direct searches, constraints 
from B(b — > S7) are very sensitive to details of the model. A small amount of flavor mixing in 
the squark sector could lead to large differences in the predictions of b — > 57 decays. For this 
reason, the constraints from B(b —>■ 57) should be interpreted with some care. 
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The large values of R in frame a) and b) in Fig. 2 encourage us to do further exploration 
in fi < case. Therefore, we plot a contours of R values in Aq - m 1( / 2 plane in Fig. 3 
for fi < case. We choose m = 300 GeV and in frame a), tan/3 = 45, while in frame 
b), tan/5 = 42. The dark-shaded regions are excluded on theoretical grounds because 
the overall theory does not lead to electroweak symmetry breaking. The slant-hatched 
region are excluded because of Z 1 is not the lightest supersymmetric particle (LSP). If 
the neutralino is the LSP, it will be stable and therefore a dark matter candidate. Within 
the slant-hatched region, charged sparticle is the LSP which disagrees with cosmologies 
models. Below the dashed line labeled "m^ = 113 GeV" are the regions is smaller 
than its experimental bound 113 GeV. Below the dotted line, branching fractions of 
B s — > are larger than 5 x 10~ 10 , but smaller than 10~ 7 . The later is the maximum 

limit the experiment, Tevatron or B factories, can explore in detecting B s — > fx + fi~ 
decay. The contours of R = B(b — > sfi + fi~)/B(b — > se + e~) are labeled by the values of 
corresponding ratio R. From frame a) we see that in the allowed region, R can be larger 
than 2. Even in frame b), R can be as large as 2. The outermost curves is for R equals 
to 1.1. Turning our attention to the sensitivity with respect to tan/5, we show in frame 
b), the same contours, but with tan/5 = 42. Although the range for tan/5 = 42 is smaller 
than that for tan/5 = 45, it has still regions where B factories may be expected to detect 
it. The ratio can be significantly larger below the R = 2 contour. 

Recently, Belle experiment has observed a signal for exclusive B decays with branching 
fraction B (B -> Kl+l~) = (0.75tgf±0.09) x 10~ 6 averaged over electron and muon chan- 
nels and B(B -> Kfi+fi') = (0.99±g;f^±g;^) x 10~ 6 in muon channel [36]. This is consistent 
with our SM value for the inclusive rate ~ 6.4 x 10~ 6 . The inclusive B decay branching 
ratio is measured by Belle group with £(5 -+X s l + r) = (6.1±lA(stat)±{i(syst))xlO- 6 
[37,38], which is consistent with the SM prediction. We also estimate that approximately 
~ 10 7 B mesons are required to reach the R = 1.2 limit. 

Due to the mf dependence of these branching fractions, one would clearly expect 
contributions to b — > st + t~ [29,39] two orders of magnitude greater than b — > sfi + fi~. 
Turning to Fig. 4, we see that this is indeed the case. In Fig. 4, we show the contours 
of B(b — > st + t~) in m — plane with A = GeV. In frame a) fj, < and tan/5 = 
45, while in frame b) /j, > and tan/5 = 53. As in Fig. 3, the dark-shaded regions 
are excluded by theoretical constraints: charge-breaking minimal or lack of electroweak 
symmetry breaking. Z\ is not the lightest supersymmetric particle (LSP) in the slant- 
hatched region. The contours are labeled by the corresponding values of b — > st + t~ 
branching fraction. The outermost contour corresponds to a branching fraction of 10~ 6 . 
The SM branching fraction is around B(b — > st + t~) ~ 3.66 x 10~ 7 . In Fig. 4, at the most 
parameter space, SUSY contributions is much larger than the SM background, reaching 
a limit of ~ 10 -4 within the allowed region. Here we remind the readers that some region 
perhaps is excluded by Bs — > fi + fi~ constraint. In ref. [9], we give out detailed calculations 
on B s —> fi + fi~. Unfortunately r identification is very difficult. The branching fraction 
about 10~ 3 for r process may be accessed at B factory. Clearly then, at least one order of 
magnitude improvement in the acceptance for b — > st + t~ is required to have a significant 
impact on this parameter space. It is interesting to note that any new physics which 
contribute significantly to b — > st + t~ will also give a signal in B s — > tt as discussed in 
Ref. [9]. 
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Another aspect of new physics contributions to the b — > sfi + fi~ signal is the forward- 
backward lepton asymmetry Afb- This is considered in ref. [25] and found to be on the 
order of a few percent hence difficult to observe at current B factory luminosities. We 
do not explore this feature in this paper, but it is potentially another channel for new 
physics detection. 

IV. SUMMARY AND CONCLUSIONS 

We have explored the rare semileptonic decay b — > se + e~, b — > sfi + fi~ and b — > 
ST + T - w ithin mSUGRA framework where some Higgs-mediated SUSY contributions to 
the branching fraction of b — > are proportional to corresponding lepton mass squared. 
When such contributions in the muon case become comparable to the SM, this opens a new 
channel for new physics detection. We explored these decays in detail and find that while 
in /j, > case, SUSY contribution is suppressed, for negative \x with large tan/5 the flavor 
violating couplings of neutral Higgs bosons to down type quarks leads to substantial SUSY 
contributions to the branching fraction of decay b — > s/i + /i~. In particular observably large 
contributions to the branching fraction of b — > sfi + fi~ are possible when tan/3 > 40. In 
particular there can be leading a large discrepancy between the b — > sfi + fi~ and b — > se + e~ 
branching ratios. The resultant deviations in R = B(b — > sfi + pL~)/ B(b — > se + e~) from 
the SM values in some regions of the mSUGRA parameter space are as shown in Fig. 3. 

In this framework the large value of m T would give a huge deviation from the SM 
prediction for B(b — > st + t~) in the region of large tan/5. We explore the branching 
fraction for b — > st + t~ decay in m — TO1/2 plane and in the region of parameter space 
considered the maximum values of B{b — > st + t~) can be up to 8 x 1CT 4 . Unfortunately, 
because of the difficulties in separating signals of r decay from the backgrounds, the r 
signal is undetectable at B factories. Thus, the most promising signal of the new physics 
considered in this paper is a deviation in the value of R from its SM value. Conversely, 
if such a deviation in R were observed it would strongly indicate new physics coupling to 
mass through the Higgs sector as in the case of mSUGRA model considered in this paper. 
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FIG. 1. One loop SUSY contributions. 
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FIG. 2. Ratio R = B(b — > s^ + fi")/B(b -► se+e") in mSUGRA model with m 1/2 = 300 GeV 
and a) R vs tan/3 with /i < 0, mo = 300 GeV, b) R vs mo with \i < 0, tan/3 = 45, c) R vs tan/3 with 
/i > 0, m = 300 GeV, d) R vs m with // > 0, tan/3 = 53. The solid line is for A Q = -300 GeV, 
and Aq = GeV for dashed line. The dotted line is the corresponding SM line. The circles 
mark the limits of the experimentally allowed regions discussed in the text. 
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m = 300 GeV, tan|3 = 45, \i < 




A (GeV) 



FIG. 3. Contours of constant ratio R = B{b — ► sfi + B{b — > se + e~) with 
|U < 0, mo = 300 GeV in — m, 1 / 2 plane with a) tan(3 = 45 b) tanfl = 42 in mSUGRA 
model. The dark-shaded region is excluded by the theoretical constraints on electroweak sym- 
metry breaking (EWSB). and the slant-hatched region is excluded for Z\ is not LSP. Below the 
dashed line is for mh < 113 GeV. Below the dotted line, the branching fractions of B s — > 
decay are larger than 5 x 10~ 10 . 
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FIG. 4. Contours of constant branching fraction for decay b — ► st + t~. The results are 
showed in mo — mi/ 2 plane with ^0 = GeV, in frame a). tan(3 = 45, fi < and in frame b). 
tan(5 = 53, \i > 0. The dark-shaded region is excluded by the theoretical constraints on EWSB. 
and the slant-hatched region is excluded for Z\ is not LSP. 
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